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ABSTRACT 


Gibbons and siamangs (Hylobatidae) are well-known for 
their rapid chromosomal evolution, which has resulted 
in high speciation rate within the family. On the other 
hand, distinct karyotypes do not prevent speciation, 
allowing interbreeding between individuals in captivity, 
and the unwanted hybrids are ethically problematic as all 


gibbon species are endangered or critically endangered. 


Thus, accurate species identification is crucial for captive 
breeding, particularly in China where studbooks are 
unavailable. Identification based on external morphology 
is difficult, especially for hybrids, because species 
are usually similar in appearance. In this study, we 
employed G-banding karyotyping and fluorescence in 
situ hybridization (FISH) as well as a PCR-based 
approach to examine karyotypic characteristics and 
identify crested gibbons of the genus Nomascus from 
zoos and nature reserves in China. We characterized 
and identified five karyotypes from 21 individuals of 
Nomascus. Using karyotypes and mitochondrial and 
nuclear genes, we identified three purebred species 
and three hybrids, including one F2 hybrid between N. 
gabriellae and N. siki. Our results also supported that 
N. leucogenys and N. siki shared the same inversion on 
chromosome 7, which resolves arguments from previous 
studies. Our results demonstrated that both karyotyping 
and DNA-based approaches were suitable for identifying 
purebred species, though neither was ideal for hybrid 
identification. The advantages and disadvantages of 
both approaches are discussed. Our results further 
highlight the importance of animal ethics and welfare, 
which are critical for endangered species in captivity. 
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INTRODUCTION 


The rates of chromosomal evolution are an order of magnitude 
higher in Mammalia than in most other classes of vertebrates 
and are strongly correlated with high speciation rates (Bush 
et al., 1977). High karyotypic diversity and rapid speciation 
are often observed in taxa characterized by small effective 
population sizes, which includes the gibbons and siamangs of 
the small ape family Hylobatidae (Primate). These animals are 
endemic to southern China, South Asia, and Southeast Asia, 
and have been assigned into four genera (Hoolock, Hylobates, 
Nomascus, and Symphalangus; Brandon-Jones et al., 2004; 
Roos & Geissmann, 2001). While siamangs (Symphalangus) 
possess a unique and large gular sac (throat pouch), gibbons 
share very similar external morphology, and were long 
acknowledged as a single genus (Hylobates). | However, 
Hoolock and Nomascus were subsequently recognized as full 
genera based on karyotypic studies, which revealed a unique 
diploid number (2n) of chromosomes among each genus: 
Hoolock (2n=38), Hylobates (2n=44), Nomascus (2n=52), 
and Symphalangus (2n=50) (Chiarelli, 1972; Prouty et al., 
1983; Van Tuinen & Ledbetter, 1983). This subdivision, as 
supported by molecular phylogenetic studies, occurred in the 
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Late Miocene approximately 7 Ma (Fan et al., 2017; Springer 
et al., 2012). Today, 20 species of gibbons and siamangs are 
recognized, though the number continues to increase due to 
newly discovered species (Fan et al., 2017). 

Speciation in this family is considerably faster than that 
observed in many other groups of mammals and is considered 
to be the result of high frequency of chromosomal rearrangement 
(Misceo et al., 2008), including compound inversion/translocation 
(Jauch et al., 1992; Koehler et al., 1995a, 1995b; Marks, 
1982; Nie et al., 2001; Stanyon & Chiarelli, 1982, 1983; 
Van Tuinen & Ledbetter, 1983; Yu et al., 1997). Previous 
comparison between the gibbon (N. leucogenys) and human 
genomes revealed 100 synteny breakpoints, which can facilitate 
chromosome breakage and rearrangement (Carbone et al., 2006). 
Gibbons exhibit highly structured society in their populations, 
which may promote chromosome rearrangement fixation, thereby 
allowing rapid speciation (Wilson et al., 1975). 

The crested gibbon genus Nomascus comprises seven 
recognized species, including N. concolor (western 
black-crested gibbon), N. gabriellae (yellow-cheeked 
gibbon), N. hainanus (Hainan gibbon), N. leucogenys 
(northern white-cheeked gibbon), N. nasutus (eastern 
black-crested gibbon), WN. siki (southern white-cheeked 
gibbon) (Brandon-Jones et al., 2004; Roos & Geissmann, 
2001), and the recently discovered N. annamensis (northern 
buffed-cheeked gibbon) (Thinh et al., 2010). These species 
are characterized by a variety of morphological, anatomical, 
karyological, and vocal features (Couturier & Lernould, 1991; 
Garza & Woodruff, 1994; Schilling, 1984; Thinh et al., 2011). 

Although all studied Nomascus species have a unique 
diploid number of 52 (Couturier et al., 1982; Dutrillaux 
et al., 1975; Van Tuinen & Ledbetter, 1983), they also 
exhibit considerable interspecific karyotypic variation. Using 
the R-banding technique, Couturier & Lernould (1991) 
distinguished the karyotypes for six individuals of N. 
leucogenys, two N. siki, three N. gabriellae, and several hybrids 
using a combination of an inversion on chromosome 7 and a 
translocation between chromosomes 1 and 22. The inversion 
on chromosome 7 was also confirmed based on the number of 
hybridization signals provided by human chromosome probes 
4 and 22. On unrearranged chromosome 7, one hybridization 
signal of probe 4 on 7q and one signal of probe 22 on 7p were 
detected. On rearranged chromosome 7, two signals of each 
probe were detected, one on each arm of the chromosomes 
(Koehler et al., 1995b). Heterozygous translocation t (1; 22) 
has also been identified in hybrids of N. leucogenys and N. 
gabriellae (Couturier et al., 1982; Turleau et al., 1983). Using 
a PCR-based approach, Carbone et al. (2009) detected 
the inversion on chromosome 7 in N. leucogenys but not in 
N. siki. These results disagreed with those of Couturier & 
Lernould (1991), and thus need to be revisited and clarified 
using additional samples. 

Currently, most gibbon species are endangered or critically 
endangered. Several species, including Hylobates lar and 
N. leucogenys, may have been extirpated from China largely 
due to poaching.  Reintroducing animals from zoos to 


their original habitat is a goal for conservation but requires 
accurate identification of purebred animals. Extreme caution 
must be undertaken in captive breeding as interspecific and 
intergeneric hybridization can occur between gibbons and 
siamangs (Hirai et al., 2007; Myers & Shafer, 1979), which 
do not resulting in reduced fertility among offspring in certain 
cases (Van Tuinen et al., 1999). Hybrid offspring are difficult to 
visually identify due to their similar external morphology, with 
identification using the karyotypic approach considered more 
reliable. Co-housing of different species has been avoided in 
many countries and animals are karyotyped for identification, 
with studbooks maintained as a database for appropriate 
management. However, this system is not yet well established 
in China and different species are usually housed together 
resulting in hybridization and unknown interbreeding history. 

In the current study, we examined 21 gibbons using a 
combination of G-banding karyotyping and fluorescence in situ 
hybridization (FISH). We also sequenced mitochondrial and 
nuclear fragments across a known inversion breakpoint on 
chromosome 7 identified in a previous study (Carbone et al., 
2009). We examined whether the combination of G-banding 
and FISH is a reliable approach for species identification. 
We evaluated the efficiency and accuracy of karyotyping- 
and PCR-based approaches for the identification of purebred 
animals and hybrids. We also re-examined whether the 
pericentric inversion on chromosome 7 is species-specific in 
N. leucogenys. 


MATERIALS AND METHODS 


Samples 

We collected peripheral blood from 20 gibbon individuals raised 
in Nanning Zoo, Ningbo Zoo, Kunming Zoo, and Huanglianshan 
National Natural Reserve, and obtained one sample (No.14 in 
Table 1) of cultured gibbon cells maintained at the Kunming 
Institute of Zoology, Chinese Academy of Sciences. Fourteen 
individuals were female and seven were male. All experimental 
procedures and animal care were performed according to the 
protocols approved by the Ethics Committee of the Kunming 
Institute of Zoology, Chinese Academy of Sciences. 


Cell culture, metaphase preparation, and G-banding 
Chromosome suspensions were obtained from lymphocyte 
cultures. Cell culture and metaphase preparations were 
performed following conventional procedures (Hungerford, 
1965). Briefly, whole blood was cultivated in the presence of 
phytohaemagglutinin in RPMI 1640 medium supplemented with 
10% fetal bovine serum. After 68—70 h of growth, colchicine 
was added to the cell cultures to a final concentration of 
0.4—0.8 ug/mL. The cell cultures were incubated for another 
2—4 h before harvest. We treated the cells with hypotonic 
solution (0.075 mol/L KCl) for 20 min, with thrice fixation in 3:1 
methanol/glacial acetic acid. Slides were prepared by applying 
10 uL of metaphase suspension onto dry and clean slides, 
which were then allowed to air dry. G-banding was performed 
following Seabright (1971). All karyotypes were analyzed after 
G-banding. 
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FISH, image capture, and processing 

We used a biotin-labeled probe of the human 22 chromosome. 
Fluorescence in situ hybridization, detection, image capture, 
and processing were carried out following Yang et al. 
(2000) and Nie et al. (2002). We detected fluorescence 
signals using a layer of Cy3-avidin (1:1 000 dilution; 
Amersham Pharmacia Biotech, USA). After detection, slides 


were mounted in Vectashield mounting medium with DAPI 
(4’6-diamidino-2-phenylindole, Vector Laboratories, USA). 
Digital images were acquired using a CytoVision system 
(Applied Imaging Co., USA) with a CCD camera mounted 
on a Zeiss microscope (Germany). We associated the 
hybridization signals with specific chromosome regions based 
on DAPI-banding patterns. 


Table 1 Samples used in this study and a summary of karyotypic characteristics including the FISH signals on the chromosome 7, 


the lengths of chromosomes 1 and 22, the results of PCR using primer sets BP and HSA, as well as the identification results 


based on each approach 

















FISH 
Nick Name House (signal on chr. 7) Chr. 1 Chr. 22 Karyotype BP HSA Cyt b genes Final identification 
(pairs) 
L h 
Su-Su Nanning Zoo 2 ong: -Short N. leucogenys N/A N/A N/A N/A 
(Norma Normal 
hgh s Long Short 
Bei-Li Nanning Zoo 2 N. leucogenys + — N. leucogenys N. leucogenys 
(Norma Normal 
3 Long Short 
Gou-Dan Nanning Zoo 2 N. leucogenys + — N. leucogenys N. leucogenys 
(Norma Normal 
L h 
San-Mei Nanning Zoo 2 ong.. Snort N. leucogenys N/A N/A N/A N/A 
(Norma Normal 
ae Long Short 
No14 Missing 2 N. leucogenys + - N. cf. leucogenys N. cf. leucogenys 
(Norma Normal 
; Long Short 
No name Kunming Zoo 2 N. leucogenys N/A N/A N/A N/A 
(Norma Normal 
L h 
NB1 Ningbo Zoo 2 ong.. Snort N. leucogenys N/A N/A N/A N/A 
(Norma Normal 
L hor 
NB2 Ningbo Zoo 2 ong; -Short N. leucogenys N/A N/A N/A N/A 
(Norma Normal 
Huanglianshan 
9 Long Short 
HLS3 Nature 2 N. leucogenys N/A N/A N/A N/A 
(Normal Normal 
Reserve 
Huanglianshan 
g Long Short 
HLS4 Nature 2 N. leucogenys N/A N/A N/A N/A 
(Norma Normal 
Reserve 
Fang-Fang Nanning Zoo 2 Short Long N. siki N/A N/A N/A N/A 
E’gui Nanning Zoo 2 Short Long N. siki + -= N. siki N. siki 
Qingguangyan Nanning Zoo 2 Short Long N. siki N/A N/A N/A N/A 
Laoer Nanning Zoo 2 Short Long N. siki N/A N/A N/A N/A 
N. siki 3 x 
Da-Shan Nanning Zoo 2 Short Long N. siki + — N. gabriellae (N. siki 3 
xN. gabriellae 2) 2 
316 Nanning Zoo 1 Short Long N. gabriellae - + N. gabriellae N. gabriellae 
Bai-Shou Nanning Zoo 1 Short Long N. gabriellae — + N. gabriellae N. gabriellae 
Jing-Jing Nanning Zoo 1 Short Long N. gabriellae — + N. gabriellae N. gabriellae 
1 short, 1 short, i N. leucogenys 
Mei-Mei Nanning Zoo 2 N. leucogenysx N. siki + - N. leucogenys A g : fae 
1 long 1 Long N. siki 
1 N. gabriellae 
Xiao-Xiao Nanning Zoo ang Short Long N. gabriellaex N. siki + + N. leucogenys g . yi i 
a half N. siki 
1 
A-Huang Nanning Zoo a Short Long N. gabriellaex N. siki N/A N/A N/A N/A 
a hal 





The karyotype and mitochondrial of Da-Shan suggest different specific affinities. +: The positive result of PCR using primer set BP or HAS. —: 


The negative result of PCR using primer set BP or HAS. N/A: Not available. 
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Amplification, sequencing, and species identification 


After karyotyping, we extracted total DNA from the cultured cells 
using a commercial DNA extraction kit (Blood & Cell Culture 
DNA Mini Kit, Qiagen, Germany). The DNA extracted from 10 
samples was considered acceptable for subsequent analyses as 
other samples were not well preserved after karyotyping. We 
amplified and sequenced cyt b genes using the primer pair 
L14724_hk3 = (5’-GGACTTATGACATGAAAAATCATCGTTG-3’) 


and H15915_hk3 (5’-GATTCCCCATTTCTGGTT TACAAGAC-3’). 


We also conducted PCR using the primers provided in Carbone 
et al. (2009) (263C9_BP_L with 263C9_BP_R (BP hereafter) and 
263C9_BP_L with 263C9_HSA22_R (HSA hereafter)), targeting 
a fragment across a breakpoint on chromosome 7. According 
to Carbone et al. (2009), HSA primers should amplify a 
fragment across a breakpoint in species without an inversion 
(e.g., Hylobates spp., N. gabriellae, and N. siki), and BP primers 
should result in positive amplification only for N. leucogenys. The 
mitochondrial and nuclear amplicons were sequenced using a 
ABI 3730 sequencer (Applied Biosystems, USA). 

We identified species based on the karyotyping and 
sequencing results using the following strategy: (1) we 
identified the species based on its karyotype; (2) we verified 
the karyotyping results using PCR with primer sets BP and 
HSA; (3) we BLAST each cyt b gene against the GenBank 
nucleotide database; (4) we estimated a gene tree to verify 
the BLAST results using our cyt b genes accompanied by 
a set of sequences representing the recognized Nomascus 
species downloaded from GenBank; and (5) we repeated 
the karyotyping and amplification/sequencing in triplicate in 
cases where the mitochondrial, nuclear sequencing, and 
karyotyping results conflicted (e.g., sample Da-Shan, see 
Results and Discussion). A final purebred or hybrid and 
species identification was determined. 

We constructed the maximum likelihood gene tree 
using RAxML v8.2.10 and the CIPRES Science Gateway 
(Stamatakis, 2014). In addition to the obtained cyt b 
sequences, 17 cyt b sequences representing seven recognized 
species were downloaded from GenBank (Supplementary 
Table S1) and aligned with our sequences using MAFFT v7.3 
(Katoh & Standley, 2014). We partitioned the alignments 
by codon positions (three partitions). We used GTR+G as 
the evolutionary model for each partition as RAxML does 
not accept models other than GTR or GTR+G. We ran each 
analysis using the rapid bootstrapping algorithm and let RAxML 
halt bootstrapping automatically. 


RESULTS 


Identified karyotypic species 


The diploid number (2n) of all analyzed gibbons was 52. 


G-banding revealed different lengths of chromosomes 1 and 
22, and FISH revealed two, three, or four fluorescence signals 
on chromosome 7 (Figure 1). Considering the G-banded 
karyotyping and FISH results, we recognized five karyotypes 
representing three karyotypic species and two hybrids, with the 
latter two distinguished by nonhomologous pairing. 
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Figure 1 Comparison of G-banded chromosomes 1, 7, and 22 
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Hybridization of HSA 22 
probe on Chr. 7 


and FISH results with human 22 chromosome-specific probe 
in different Nomascus species and their hybrids 

A: Nomascus leucogenys; B: Nomascus siki; C: Nomascus gabriellae; D: 
hybrid of N. leucogenysxN. siki; E: hybrid of N. gabriellaex N. siki. 


Nomascus leucogenys (Figure 2) 

The G-banded karyotypes of 10 individuals (7°, 33) were similar 
to the G-banded (Van Tuinen & Ledbetter, 1983) and R-banded 
karyotypes for N. leucogenys (Couturier & Lernould, 1991). 
Chromosomes 1 and 22 were normal. Two pairs of FISH signals 
were found on chromosome 7, supporting a pericentric inversion 
event (Figure 1A). 
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Figure 2 G-banded karyotype of a northern white-cheeked 
gibbon (N. leucogenys) 


Nomascus siki (Figure 3) 


The karyotypes of five samples (3, 28) were identical to 
those of N. siki reported by Couturier & Lernould (1991). 
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Compared with N. leucogenys, a balanced translocation 
between chromosome pairs 1 and 22 was detected, resulting 
in one pair of shortened chromosome 1 and one pair of derived 
chromosome 22 (t (1; 22)). Similar to that observed in N. 
leucogenys, FISH demonstrated a pericentric inversion event 
on chromosome 7 (Figure 1B). The sample “Da-Shan” had a 
typical siki-karyotype. 
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Figure 3 G-banded karyotype of a southern white-cheeked 
gibbon (N. siki) 


Nomascus gabriellae (Figure 4) 

The G-banded karyotypes of three animals (19, 24) were the 
same as that of N. gabriellae (Couturier & Lernould, 1991). 
This karyotype was similar to N. siki in the presence of t (1; 22). 
Only one pair of FISH signals were observed on chromosome 
7, thus differing from N. siki and N. leucogenys (Figure 1C). 
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Figure 4 G-banded karyotype of a yellow-cheeked gibbon (N. 
gabriellae) 


Nomascus leucogenysx Nomascus siki hybrid (Figure 5) 
The G-banded karyotype of one specimen was unique, 
characterized by a heterozygous translocation. One 


chromosome 1 and 22 were normal, similar to that of N. 


leucogenys, and the other chromosome 1 and 22 were similar 
to that of N. siki, indicating one translocation t (1; 22). Two pairs 
of FISH signals were observed on chromosome 7, indicating 
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a pericentric inversion event (Figure 1D). This specimen was 
identified as a hybrid of N. leucogenysx N. siki. 
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Figure 5 G-banded karyotype of a N. leucogenysxN. siki 
hybrid 


Nomascus gabriellaex Nomascus siki hybrid (Figure 6) 
The distinct G-banded karyotype of two individuals was similar 
to that of N. siki and N. gabriellae in the presence of 
t (1; 22). FISH detected three fluorescence signals on 
chromosome 7, indicating heterozygous pericentric inversion 
(Figure 1E). These two individuals were identified as hybrids of 
N. gabriellaex N. siki. 
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Figure 6 G-banded karyotype of a N. leucogenysxN. siki 
hybrid 


Mitochondrial and nuclear sequences 

We determined the complete cyt b sequences for 10 individuals 
representing the three purebred species (n=7) and two 
hybrids (n=3) identified in the karyotypic analyses. All new 
sequences are available in GenBank under accession Nos. 
MH188408—MH188428 (Supplementary Table S1). Five out 
of the seven purebred animals were supported by BLAST and 
phylogenetic analyses using the cyt b gene and the other two 
were not (Figure 7). One sample (Da-Shan) was identified 
as N. siki based on its karyotype, but its mitochondrial gene 
was typical of N. gabriellae, which also supported by our 
phylogenetic analyses (bootstrap value (BS)=98). Sample 


14 was identified as N. leucogenys based on its karyotype 
but could not be identified unambiguously as it was equally 
similar to N. leucogenys and N. siki. It was closely related 
to (BS=84) and formed a sister lineage of these two species 
on the phylogenetic tree (BS=52). The cyt b genes of the 
two hybrids identified in the karyotypic analyses (Xiao-Xiao 
and Mei-Mei) corresponded to one of their parental species 
(BS>88). On our cyt b genetic tree (Figure 7), sample 14 was 
a close relative to (BS=84) and formed a sister lineage with N. 
siki+N. leucogenys (BS=52). 
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Figure 7 Maximum-likelihood gene tree estimated using the 
cyt b genes of Nomascus species 
Branch lengths correspond to substitution rates. Numbers on branches are 
bootstrap values. Sequences downloaded from GenBank are in gray and 
GenBank accession Nos. are shown. 


The PCR results using the HSA and BP primer sets were 
concordant with the karyotypic results. The PCR results using the 
BP primers (confirming the inversion) were positive for purebred 
N. leucogenys, N. siki, and the N. leucogenysxN. siki hybrid 
(Mei-Mei) and were negative when using the HSA primers for 
these samples (Table 1). The PCR results using the BP primers 
were negative for purebred N. gabriellae but positive when using 
the HSA primers. For the hybrid identified as N. gabriellaex N. 
siki (Xiao-Xiao), PCR and sequencing were successful using 
both BP and HSA primers. For Da-Shan, the PCR results were 
positive using BP primers and negative using HSA primers, which 
were congruent with its karyotype (siki type) but disagreed with 


its mitochondrial gene results. 

Because the karyotype and PCR results were congruent, we 
assigned seven samples as purebred species and identified 
another two as hybrids. One sample (14) was identified 
as N. cf. leucogenys as it did not cluster with the other 
N. leucogenys sequences downloaded from GenBank. We 
considered Da-Shan as an F2 hybrid of N. gabriellae and N. 
siki due to the inconsistent patterns revealed by karyotyping 
and PCR using the HSA and BP primers (siki), and the 
mitochondrial phylogeny (gabriellae) (Table 1). 


DISCUSSION 


We used karyotypic- and PCR-based approaches to examine 
the affinities of Nomascus species characterized by high 
karyotypic diversity. We determined the specific karyotypes 
of each species as well as the existence of an inversion on 
chromosome 7 in N. siki and N. leucogenys. We also revealed 
the definitive existence of hybrids in zoos in China, which calls 
attention to the animal ethics and welfare issues related to 
endangered species breeding in captivity. Using the Nomascus 
species as our focal taxa, we evaluated the accuracy and 
efficiency of both approaches in identifying purebred and hybrid 
species. 

Our results supported that N. gabriellae, N. leucogenys, and 
N. siki are distinguishable based on translocation t (1; 22) 
and inversion on chromosome 7 (Figures 1—5), even though N. 
leucogenys and N. siki diverged from each other only recently. 
Couturier & Lernould (1991) revealed that N. hainanus does not 
have an inversion on chromosome 7 or translocation t (1; 22), 
and therefore has a distinct karyotype. It would be interesting 
to examine the G-banded karyotypes and FISH results of 
N. concolor as well as N. annamensis, which is sister to N. 
gabriellae and recognized only recently (Thinh et al., 2010). 
Due to interspecific karyotypic variation, hybrids can be easily 
distinguished based on heterozygous R-, G-banded karyotypes 
and FISH results (Couturier & Lernould, 1991; this study). In 
wild populations, hybridization of Hylobates species such as H. 
larx H. pileatus, H. larx H. agilis, and H. muellerix H. agilis/H. 
albibarbis has been observed (Brockelman & Gittins, 1984). 
Hybridization between Nomascus species in the wild has not 
been reported, which is likely due to their allopatric distribution. 
However, we observed heterozygous karyotypes in captive 
Nomascus animals, congruent with previous study (Couturier 
& Lernould, 1991). Our results support the possibility of 
nonhomologous pairing during meiophase (Figures 5,6), which 
is an interesting characteristic in Hylobatidae. Because 
gibbons have very strong social structures and usually live 
in small family groups, such nonhomologous pairing may be 
easily fixed during evolutionary histories and further promote 
diversification and speciation, resulting in high species diversity 
within the family. 

Our FISH and DNA analyses supported the existence of 
an inversion on chromosome 7 in both N. leucogenys and N. 
siki (Figures 1—3), consistent with the findings of Couturier 
& Lernould (1991), though not supported by Carbone et al. 
(2009). The negative results obtained in the latter study may 
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be due to the hair sample used to represent N. siki, which is 
known for low DNA quantity as well as the existence of PCR 
inhibitors. 

The phylogenetic position of sample 14 was an interesting 
finding. Its karyotype was identical to the other N. leucogenys 
samples (Table 1), which was not supported by the 
mitochondrial gene tree (Figure 7). We repeated amplification 
and sequencing three times for this sample. The sequences 
were identical, indicating no cross-sample contamination, 
and no premature codon was observed, indicating it was 
not a pseudogene. There are two hypothetical scenarios 
that may explain these results. One is incomplete lineage 
sorting between N. leucogenys and N. siki, resulting in 
non-monophyletic relationships. This situation has never been 
observed because the effective population sizes of gibbon 
species are usually small. Alternatively, this sample might 
represent a distinct and unknown taxon, which is a close 
relative but not identical to N. leucogenys or N. siki. The 
sample was obtained from the Kunming Institute of Zoology, 
Chinese Academy of Sciences, in 1993 and the information on 
the animal was not recorded in detail. Unfortunately, neither 
hypothesis could be supported or rejected in the current study. 

Similar to previous studies, our study supported that 
hybridization has occurred in captivity. However, fertility may 
be less impacted as identification of Da-Shan showed it to be 
a likely F2 hybrid of N. gabriellae and N. siki, despite there 
being no available studbook. This is the first record of a 
cross-back F2 hybrid in Nomascus as most other countries 
have well-established systems in place to prevent interbreeding 
of gibbons in captivity, with many animals also previously 
karyotyped. Continuous interbreeding and production of fertile 
offspring may spread across zoos in China because there are 
only limited numbers of gibbons in zoos, and correct species 
identification of their offspring will become far more difficult 
after several generations and recombination. Our findings call 
for the introduction of a system to prevent gibbon interbreeding 
in captivity and for better welfare and awareness of these 
animals. 

The R-banding technique can distinguish different 
karyotypes of Nomascus species (Couturier & Lernould, 
1991), but is both difficult and time-consuming. Herein, we 
demonstrated that G-banding in combination with FISH is 


a reliable approach for identifying karyotypic characteristics. 


This approach was also appropriate for species identification, 
though it was limited in identifying cross-back hybrid F2 


individuals and did not recognize Da-Shan as a hybrid (Table 1). 


As per Carbone et al. (2009), we agree that PCR and 
sequencing can be applicable for examining chromosome 
inversion and species identification. This approach does not 
require high quality cultured cells or karyotyping techniques 
and is applicable for DNA samples with low quality and/or 
low yield, which certainly include, but are not limited to saliva, 
urine, hair, and feces. The known limitation is that accurate 
identification of hybrids relies on a finely established system 
with known karyotypes, breakpoints, and primers. In our case, 
it easily distinguished N. leucogenys/N. siki from N. gabriellae, 
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but could not easily distinguish N. leucogenys/N. siki from the 
hybrid of N. leucogenysx N. siki. 
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